Longitudinal spin decoherence in spin diffusion in semiconductors 
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We set up a set of many-body kinetic Bloch equations with spacial inhomogeneity. We reexamine 
the widely adopted quasi-independent electron model and show the inadequacy of this model in 
studying the spin transport. We further point out a new decoherence effect based on interference 
along the direction of diffusion in spin transport due to the so called inhomogeneous broadening 
effect in the Bloch equations. We show that this inhomogeneous broadening can cause decoherence 
alone, even in the absence of the scattering and that the resulting decoherence is more important 
than the dephasing effect due to both the D'yakonov-Perel' (DP) term and the scattering. 



PACS: 72.25.Dc; 72.25.Rb 

Study of spintronics has attracted tremendous atten- 
tion in recent years, both in theoretical and experimen- 
tal circles,El thanks to the discovery of the long-lived 
(sometimes > 100 ns) ciiherjanL electron spin states in 
n-typed semiconductors.HuQQ'Q'tj Possible applications of 
spintronics include qubits for quantum computers, quan- 
tum memory devices, spin transistors, and spin valves 
etc. The last two applications involve transporting spin 
polarized electrons from a place to another by means 
of an electrical or diffusive current. Therefore, it is of 
great importance to study the spin transport. Apart 
from the great number of works on spin injection, there 
are only a few experimental reports, sm coherent spin 
transport over macroscopic distanccotrEl On theoretical 
aspect, most works are based on a quasi-independent 
electronu-jmad^ and .focused on the diffusive transport 
regime ,^00'E3 , EJ'Ej where equations for spin polarized 
currents can be set up and the longitudinal spin dephas- 
ing, generally referred to as spin diffusion length can be 
achieved. In these theories, the mechanism for the spin 
relaxation is assumed to be due to the spin-flip scatter- 
ing. In the absence of the scattering, the spin polariza- 
tion will not decay in a nonmagnetic sample. In Ref. |l6| , 
Takahashi et al. calculated the scattering induced spin 
relaxation time associated with the spin diffusion starting 
from the many body kinetic equations. 

Of particular interest to the spin transport theory 
in semiconductors has been the question as to whether 
the quasi-independent electron model can adequately ac- 



count for the experimental results or whether many-body 
processes are important. Flatte et al. have concluded 
that an independent electron approach is quite capable 
of explaining-measurements of spin lifetimes in the diffu- 
sive regimeO In this paper, we reexamine this issue from 
a full many-body transport theory and show the inade- 
quacy of the independent electron model in describing 
the spin transport. We also propose a mechanism that 
may cause strong longitudinal spin decoherence in addi- 
tion to the spin dephasing due to scattering. The new 
mechanism is based on the interference effect due to the 
wavevector dependence of the spin densities along the 
spacial gradients in the spin diffusion. This wavevector 
dependence can be considered as some sort of "inhomo- 
geneous broadening" , which can cause spin decay alone, 
even in the absence of scattering. 

Recently, we have presented a many-body kinetic the- 
ory to describe the spin precession and dephaajjp^-ia. in- 
sulating samples as well as n-doped samplesE§E3cI] In 
this paper we extend this theory to the spacial inho- 
mogeneous regime and obtain the many-body transport 
equations necessary to investigate the spin diffusion in n- 
doped GaAs. Here, we only focus on the spin transport 
inside the semiconductor and avoid the problem of spin 
injection at the boundary. Based on the two-spin-band 
modelij in the conduction bands, we construct the semi- 
conductor Bloch equations by using the nonequilibrium 
Green function method with gradient expansion as well 
as the generalized Kadanoff-Baym AnsatzEll as follows: 
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Here p(R>, k, t) represents a single particle density ma- 
trix. The diagonal elements describe the electron distri- 
bution functions p aa (R, k, t) = f a (R, k, t) of wave vector 
k and spin <r(= ±1/2) at position R and time t. The off- 
diagonal elements /v_ CT (R, k, t) describe the inter-spin- 
band polarization components (coherences) for the spin 



coherence. The quasi-particle energy s aa i (R, k, t) , in the 
presence of a moderate magnetic field B and with the 
DP mechanismtJ included, can be written as 

£ CTCT / (R, k, t) = £fe<W + [gpB'B + h(k)] • 



■eV>(R,t) + 5w(R,M) 



(2) 
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Here Sk — k 2 /2m* is the energy spectrum with m* 
denoting electron effective mass, — e is the electron 
charge and a are the Pauli matrices and h(k) origi- 
nate from the DP mechanism which contains both the 
DresselhauscJ and the Rashba terms .o In this paper, 
we only consider the-iirst one. For [001] quantum well, 
it can be written asc3 h x (k) — jk x (k 2 — K 2 ,), h y (\c) = 
jkyfaz — k 2 ), with k 2 denoting the average of the opera- 
tor — (d/dz) 2 over the electronic state of the lowest sub- 
band. 7 = (4/3)(™7m CT )(l/ V2m* 3 £ g )(77/Vl - v/3) 
and 77 = A/ (E g + A). Here E g denotes the band gap, 
A represents the spin-orbit splitting of the valence band, 
and m cv is a constant close in magnitude to the free elec- 
tron mass m E3 The electric potential t/j(R.,t) satisfies 
the Poisson equation 



V^R, t) = -e [n(R, t) - n (R)] /e 



(3) 



where n(R,t) = J2ek /<r(R; k, t) is the electron den- 
sity at position R and time t, and no(R) is the back- 
ground positive electric charge density. £ CTCT / (R, k, i) = 
— Vqp a(7 ' (R, k — q, t) is the Hartree-Fock self-energy, 
with Vq denoting the Coulomb matrix element. In 2D 
case, Vq is given by 



Aire 2 



~ e (q 2 +q 2 + n 2 ) 



\HiQz)\ 



(4) 



in which n = 2e 2 m*/eo fa(K = 0) i s the inverse 
screening length, with £0 being the static dielectric con- 
stant. The form factor \I(iq z )\ 2 = n 2 sin 2 y/[y 2 (y 2 — tt 2 ) 2 ] 
with y = q z a/2. It is noted that when one takes only 
the diagonal elements p aa of Eq. ([!]) and neglects all 
off-diagonal ones p a -a, the first three terms on the left 
hand side of the equation correspond to the drift terms in 
the classical Boltzmann equation, modified with the DP 
terms and self energy from the Coulomb Hartree term. 

dP ^ at'*'*' 1 lc an< i 9p ^t k ' t " > Is m the Bloch equations ([l]) are 
the coherent and scattering terms respectively, with the 
symbols | c and | s standing for "coherent" and "scatter- 
ing" . T-ha,components of the coherent terms can be writ- 
ten asfeilEl 



dt 

dPa-a 



dt 



-2Im[e cr _ -/?_ CTff ], (5) 

i [e aa - £-<t-<t\ Pa- a + i£ a -a [f-a ~ fa] ■ (6) 

) 



While the scattering terms a p(R^Mj are given in detail 

in Eqs. (5) and (7) of Ref. ||. 

The Bloch equations (|l|) can be reduced to their coun- 
terpart in the independent electron approach as follows. 
The DP term forms an effective magnetic field. It can 
flip the spin-up electrons to the spin-down ones, and vice 
versa. The DP term combines with the. scattering will re- 
sult in a longitudinal spin dephasing E3' E3l3 By applying 
the relaxation time approximation to describe this de- 
phasing and discarding the spin coherences ^-^(R, k, t) 



as well as the DP term (to avoid double counting) and 
carrying out the summation over k, one obtains the the 
continuity equation for electrons of spin a 

dn a (R,t) 1 ry(R,t)-n (R, t) 
— -V R -J CT (R,t) = , 

(7) 

in which no(R, t) = [n a (R, t) + n_ ff (R, t)]/2 is the total 
electron number at R. J<x(R, t) = ^ k (— e)v ff k/cr(R, k, t) 
is the electric current of spin a. The spin dependent ve- 
locity is Vg-k = Vi^cro^R, k, t) where e^^R, k, t) is given 
by Eq. (Q) but without the DP term h(k). By applying 
the relaxation time approximation to describe the mo- 
mentum scattering and keeping terms of the the lowest 
order (ie., neglecting terms containing p a ~a) and carry- 
ing out the summation over k, one obtains the expression 
for the current in the steady state: 

J CT (R, t) = n CT (R, t)epE(R, t) + eL>V R n ff (R, t) . (8) 

Here p and D represent the electron mobility and dif- 
fusion constant respectively. Equations (0) and (||) 
are the diffusion, aquations in the independent electron 
approachMBEyia 

One can see from the derivation of above diffusion 
equations that, by summing over k, the k dependence 
of the coefficients of Vr/o(R, k, t) in the Bloch equation 
(Q) is removed. This will not cause any problem when 
there is no spin precession. However, when the electron 
spin precesses along with the diffusion, e.g. in the pres- 
ence of a magnetic field or of an effective one (i.e. the DP 
term), this kind of k dependence may cause additional 
decoherence. 

To reveal this effect , we apply the above kinetic equa- 
tion to study the stationary state in the plane of an n- 
doped GaAs quantum well (QW) , with its growth direc- 
tion along the z-axis. The width of the QW is assumed 
to be small enough so that only the lowest subband is 
important. We assume one side of the sample (x = 0) is 
connected with an Ohmic contact which gives constant 
spin polarized injection. In this study, we assume the 
electric field E = 0. The diffusion is along the x direc- 
tion. The electron distribution functions at the interface 
are assumed to be the Fermi distributions 

U(0,k,t) = /°(k) = {exp[( £fc - p a )/T] + l}- 1 , (9) 

with T being the temperature and p a representing the 
electron chemical potential of spin a. The spin coher- 
ence at the interface is assumed to be zero 



p^ a (0,k,t) = 



(10) 



It is understood that the boundary condition here is an 
approximation to describe the distributions just after the 
injection of the spin polarization from the Ohmic con- 
tact. There is no net charge injection into the QW and 
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the well is kept charge neutral everywhere. Actually, this 
boundary condition does not necessarily come from the 
injection at the interface. It can also be produced in 
the center of semiconductors by a circularly polarized cw 
laser. 

We first consider a much simplified case by neglecting 
the DP terms h(k), the self energies as well as the scat- 
tering terms in the Bloch equations (Q). The simplified 
equations are therefore as follows 

k 

-^-d x fa(x,k) - gn B Blm[p- a . a (x,k)] = 0, (11) 

171 

%d x p CT _ CT (x,k)-i^^A/ CT (z,k)=0. (12) 
m* 2 

Here we take the magnetic field B along the x-axis. 
Af a (x,k) — f a (x,k) — /_ CT (a;,k). The solution for this 
simplified equations with the boundary conditions (||) 
and @ can be written out directly 

A/ CT (x i k)=A/ (k)cos gMB f m * 3; , (13) 

i i ^ i a fOn \ ■ 9PBBm*x 
p a -a{x, k) = -A/ U (k) sin . (14) 

Equations ( |l3| ) and (|lj) clearly show the effect of the 
fc-dependence to the spin precession along the diffusion 
direction. For each fixed k x , the spin precesses along 
the diffusion direction with fixed period without any de- 
cay. Nevertheless, for different k x the period is differ- 
ent. The total difference of the electron densities with 
different spin is the summation over all wavenumbers 
AiV = ^ k Af a (x, k). It is noted that the phase at the 
contact x = for different k x is all the same. How- 
ever, the speed of the phase of spin precession is different 



for different k x . Consequently, when x is large enough, 
spins with different phases may cancel each other. This 
can further be seen from Fig. 1 where the electron den- 
sities N a = ^ k ftr(x, k) for up and down spin are plot- 
ted as functions of position x. The boundary electron 
densities at x = are N 1 / 2 (0) = 2.05 x 10 11 cm~ 2 and 
N-i/2(0) = 1.95 x 10 11 cm" 2 . We take B = 1 T and 
T = 200 K. In order to show the transverse spin dephas- 
ing, we plot in the same figure the incoherently summed 
spin coherence p{t) = J2w \p^-^( x ^)\- It is understood 
that both the true dissipation and the interference among 
the k states may contribute to the decay. The decay 
due to interference is caused by the different precessing 
rates of electrons with different wavevectors. For finite 
system, this-Leads to reversible loss of coherence among 
electrons EZEa We refer to this kind of loss of coherence 
as decoherence. Whereas for the true dissipation, the 
coherence of the electrons is lost irreversibly.EjOEa The 
irreversible loss of coherence is termed dephasing in this 
paper. The incoherent summation is therefore used to 
isolate the ip/azersible decay from the decay caused by 
interference.E3'E3 From the figure, one can see clearly the 
longitudinal decoherence caused by the interference ef- 
fect. It is also noted from the figure that p does not 
decay with the distance. This is consistent with the fact 
that there is no scattering in Eqs. ([D]) and (|l4|) and the 
decay comes only from the interference effect. 

Facilitated with the above understanding, we turn 
to the spin diffusion problem with the DP terms, self- 
energies and scattering included. We take B = E = 0. 
By substituting the quasi-particle energy e Ga < (R, k, t) 
[Eq. ([|)] into the Bloch equations ([!]), the first three 
terms in Eqs. (Q) can be written as 



dtPaa' (R, k, t) + ed x tp(R, t)d ksr p aa > (R, k, t) - i ^ [^S,^ (R, k, tjd^p^a' (R, k, t) + d ki p arJl (R, k, t)d x Ys aia , (R, k, t) 



k 1 r 

•— docPaa'CR, M) + - d x (h x (k) - iahy(k))d x p-aar(R,k, t) +d kx {h x (k) + io'h y (\t))d x p a - a > (R, k, t) 
o [ 5fe - ( R ' k ' ( R ' k ' *) + dxp aa i (R, k, t)d kx £ CTlCT ' (R, k, t) 



(15) 



It is therefore noted that the corresponding coefficients 
of d x p a(7 i , dxP-cru' and d x p a - a > in the Bloch eqautions 
are 

% + iaU£-(R,M) + iW(R,M)] , (16) 

m* 2 

ift.{[ft x (k)-»CT/i 1 ,(k)]/2 + E D ._ .(R,k,t)} , (17) 

i^{[/i a (k)+iaX(k)]/2 + S_ ffV ,(R,k,t)} , (18) 

respectively. They are all fc-dependent. Hence, similar 
to the simplified model, the interference effect is also im- 
portant in the full kinetic equation. The kinetic equa- 



tions (0) and the Poisson equation (||) , together with the 
boundary conditions (||) and ([lo]) can be solved numer- 
ically iii|-an_iterative manner to achieve the stationary 
solution.E3E!l The numerical results for a typical QW 
with width a = 7.5 am, boundary spin polarization 
JV 1/2 (0) = 2.05 x 10 11 cm' 2 and 2V_ 1/2 (0) = 1.95 x 10 11 
cm~ 2 at temperature T = 200K are plotted in Fig. 2. In 
this computation, we only take into account the scatter- 
ing due to longitude optical (LO) phonon. It can be seen 
from the figure that the surplus of the spin up electrons 
decreases rapidly along the diffusion direction, similar to 
the simplified model shown above. 
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FIG. 1. Electron densities of up spin and down spin (solid 
curves) and incoherently summed spin coherence p (dashed 
curve) versus the diffusion length x. B — 1 T. Note the scale 
of the spin coherence is on the right side of the figure. 
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FIG. 2. Electron densities of up spin and down spin and 
the incoherently summed spin coherence versus the diffusion 
length x. Solid curves and dashed curve: N a and p from the 
full Bloch equations; Dash-dotted curves and dotted curve: 
N a and p from the equations without the interference effect. 
Note the scale of the spin coherence is on the right side of the 
figure. 

The fast decay above is understood mainly generated 
by the decoherence from the interference effect due to 
the inhomogeneous broadening. Other dephasing effects 
such as those caused by the DP terms in Eqs. (|J) and (||) 
as well as the spin conserving LO phonon scattering also 



contribute to the decay. Besides, we pointed out that the 
inhomogeneous broadening effect combined with spip* 
conserving scattering can also cause spin dephasing. t2l 
Therefore, the above mentioned inhomogeneous broad- 
ening may also cause spin dephasing in the presence of 
LO phonon scattering. To compare the decoherence due 
to interference and the dephasing due to the DP term 
together with the scattering, we remove the interference 
effect in the transport equations by replacing k in the 
coefficients [Eqs. (fjj§)-([L8|)] with k = kp. Here &f rep- 
resents the Fermi wavevector. Therefore, if there is any 
decay of spin polarization along the diffusion direction, 
it comes from the spin dephasing. The numerical result 
is plotted in Fig. 2. It is shown clearly that the decay of 
spin polarization due to the dephasing effect alone (dash- 
dotted curves) is much slower than that due to the de- 
coherence (interference) effect. In the figure we also plot 
the corresponding incoherently summed spin coherences 
p. One can see from the figure that both coherences p 
decay slowly and their decay rates are comparable when 
x > 1 /xm. This further justifies what mentioned above 
that the fast decay of the spin polarization is mainly due 
to the interference effect. 

In conclusion, we have set up many-body kinetic Bloch 
equations with spacial inhomogeneity. We reexamined 
the wildly adopted quasi-independent electron model and 
pointed out an important many-body spin decoherence 
effect which is missing in the single electron model. The 
new decoherence effect is based on an interference ef- 
fect along the diffusion direction in spin transport due 
to the so called inhomogeneous broadening effect. We 
have shown that this inhomogeneous broadening effect 
can alone cause spin decoherence, even without the scat- 
tering and that the resulting decoherence is more impor- 
tant than the dephasing effect due to both the DP term 
and the scattering. Our study shows the inadequacy of 
the quasi-independent electron model. Therefore, it is 
important to use the full many-body theory to study the 
spin transport. 
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